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Magnetic damping of the spin, the decay rate from the initial spin state to the final state, can be controlled by 
the spin transfer torque. Such an active control of damping has given birth to novel phenomena like the 
current- driven magnetization reversal and the steady spin precession. The spintronic devices based on such 
phenomena generally consist of two separate spin layers, i.e., free and pinned layers. Here we report that the 
dipolar coupling between the two layers, which has been considered to give only marginal effects on the 
current driven spin dynamics, actually has a serious impact on it. The damping of the coupled spin system 
was greatly enhanced at a specific field, which could not be understood if the spin dynamics in each layer was 
considered separately. Our results give a way to control the magnetic damping of the dipolar coupled spin 
system through the external magnetic field. 

Spins misaligned to the magnetic field relax into the field direction by transferring their energy to the 
environment. The relaxation rate is governed by the magnetic damping constant (ot) which reveals the 
interaction between the spins and the environment around the spins. Though ot is a material- specific value, 
it can be actively controlled by the spin transfer torque^'^. The active reduction of a gives such interesting 
phenomena as magnetization reversaP"^ and steady spin precession^"^^, which have given birth to the novel 
devices like the spin transfer torque magneto resistive random access memory (STT-MRAM)^^'^^ and the spin 
torque nano- oscillator (STNO)^^. In general, these spintronic devices consist of two main spin layers called as free 
and pinned layer. Although each layer is coupled to the other by the dipolar interaction and the spin transfer 
torque, the spin dynamics driven by their coupling has not been considered seriously. Only recently, the spin 
dynamics coupled by the mutual spin transfer torque has been studied by micromagnetic calculation^^'^^. But the 
dipolar interaction was assumed to give only marginal effects because the dipolar field from the pinned layer is 
generally compensated by the antiferromagnetically coupled additional layer, i.e. synthetic antiferromagnet 
(SAF)^^"^^. Here in this report, however, we show that the dipolar coupling actually has a serious impact on 
the current driven spin dynamics. We observed that the current driven spin oscillation disappears at a specific 
magnetic field. The micromagnetic calculation and the numerical estimation of the eigenmodes on the model 
systems with different dipolar interactions show that the observed behavior originates from the dipolar coupling. 
As will be described later, the current driven magnetic damping of the dipolar- coupled spin system is closely 
related to the transfer of energy between the two layers. By changing the energy transfer rate with the external 
magnetic field, one can adjust the damping of the coupled spin system. 

Results 

Anomalous discontinuity in spin oscillation mode observed in magnetic tunnel junction. The sample for our 
study is a typical magneto resistive tunneling junction (MTJ), which consists of the bottom layer [Si/Si02/ 
TiN(60)], SAF [PtMn(15)/Co9oFeio(1.5)/Ru(0.8)/Co4oFe4oB2o(1.5)], tunneling barrier [Mg(0.3)/MgO(0.53)], 
magnetic free layer [Co2oFe6oB2o(2.0)], and capping layer[Ta(5)/Ti(15)/Ru(15)] as shown in Fig. 1. Here the 
numbers in the parenthesis represent the nominal thickness in the unit of nanometer. In order to study the 
coupled motion between the free layer (FL) and the top pinned layer of the SAF (TPL), the thickness of the pinned 
layer was thin enough to be excited by the spin transfer torque. In previous studies that focus only on the dynamics 
of the FL spin, the thickness of the TPL was much thicker^°'^^ or wider^^ than the FL, inhibiting the excitation of the 
TPL. The MTJ cell has a circular shape with a diameter of 90 nm, minimizing the in-plane shape anisotropy. The 
current driven spin precession yields a resistance oscillation in the circular MTJ cell, which has been measured 
under different external fields {fioH)- The polarity of was defined as positive when the magnetization of the FL 
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Figure 1 | Schematic diagram of the MTJ structure and measurement setup. The three ferromagnetic layers are magnetically coupled. The free layer 
(Co2oFe6oB2o) is dipolar-coupled with the top pinned layer (Co4oFe4oB2o). The top pinned layer is exchange biased antiferromagnetically with the bottom 
pinned layer (CogoFcio) through the Ru layer, forming a synthetic antiferromagnet. For the study of coupled spin motion, the thickness of the pinned layer 
was thin enough to be excited by the spin transfer torque. The current driven spin oscillation yields the resistance change in the MTJ, the microwave by 
which is measured through the spectrum analyzer (S.A.). 



is antiparallel with the TPL. The current bias polarity was defined as 
positive when the electron moves from the TPL to the FL. 

Figure 2a shows the resistance of the MTJ cell as a function of /HqH 
in the current density / = 7.8xl0^ A/ cm^. Under the negative (pos- 
itive) field, the spins in the FL are parallel (antiparallel) to the spins in 
the TPL, yielding the lower (higher) magnetoresistance level. The 
spins in the TPL begin to be tilted from the spin flop-field jHoHsf = 
90 mT, where the external field becomes comparable to the exchange 
field between the TPL and bottom pinned layer (BPL) of the SAF. 
The alignments of the spins in the three layers are schematically 
shown in Fig. 2a as arrow marks for the comparison with the resist- 
ance curve. 

The color- coded power spectrums of the MTJ cell in the current 
biases / = +7.8 x 10^ A/cm^ and -7.8 x 10^ A/cm^ are shown in 
Fig. 2b and 2c, respectively. In the field range 0 < /HqH < jHoHsfi one 
can observe distinctive dependence of the oscillation mode on the 
current bias polarity. In positive (negative) current bias, the spins in 
the FL (TPL) are dominantly excited, yielding a blue (red) frequency 
shift with the increasing field. Above the spin-flop transition /HqH > 
iLioHsfi both the free and pinned layer show blue frequency shifts. 
These general features agree with the previous reports on the eigen- 
modes of the FL and TPL in the MTJ celP'-'^ 

There are several interesting points in Fig. 2b and 2c. A discon- 
tinuity is clearly observed in the oscillation mode. In Fig. 2b, the free 
layer oscillation observed at low fields disappears at around /HqH = 
40 mT, and reappears at /HqH > 60 mT. One should note that the 
point by point discontinuity in the field range /HqH < 30 mT is just 
due to the large interval of the sampling fields in our experiment. 
This is different from the real discontinuity of the oscillation mode. 

One can also find a gap in the oscillation frequency at around the 
disconnected point. The extension of the lower field mode by the 
Kittel formula^^ does not matches well with the higher field mode. 
The oscillation mode of the TPL in Fig. 2c also shows similar dis- 
continuity and frequency gap, which implies that the observed fea- 
ture is independent of the current polarity. The explanation of the 
observed anomalous breakdown of the oscillation mode is the main 
subject of this article. 

Micromagnetic calculation : effect of dipolar coupling on spin 
oscillation. The magnetization dynamics of the coupled spin layers 



have been calculated using the Landau-Lifshits-Gilbert-Slonczewski- 
Slonczewski (LLGS) equation given by^'^^'^^ 

dM a dM aj , n . n 

= —yiiMxH-\ Mx— y — Mx(Mxp) (1) 

dt M dt ^ M V i^; V ; 

where y, a, Up and p represents the gyromagnetic ratio, the intrinsic 
damping constant, the amplitude of the spin-torque in the unit of 
magnetic field, a unit vector parallel to the electron spin polarization. 
In order to study the effect of the dipolar coupling between the layers, 
the micromagnetic calculations have been carried out on several 
different separations (s) between the TPL and FL from 0.8 nm to 
50 nm. Although electrons cannot tunnel through the oxide barriers 
in case of large s value in reality, we assumed the same tunneling 
current could pass through the MTJ in our calculations. 

The color-coded microwave power calculated on the model sys- 
tem with s = 0.8 nm, 10 nm, 50 nm under the positive and negative 
current biases /= +8.0 x 10^ A/cm^ are depicted in Fig. 3a-3c and 
in Fig. 3d-3f, respectively. The case of s = 0.8 nm shows nearly the 
same oscillation feature with the experiment. There is a clear discon- 
tinuity and frequency gap in the oscillation mode at around 45 mT. 
As s increases, i.e., the dipolar interaction between the FL and TPL 
becomes weaker, the gap between the higher frequency mode and 
lower frequency mode decreases. At the separation of 50 nm, the 
oscillation mode is nearly connected with a negligible gap like the 
oscillation mode of a 'non-coupled' single spin system. This reveals 
that the observed breakdown and frequency gap in the oscillation 
mode at around 45 mT comes from the dipolar coupling between the 
FL and TPL. 

Normal modes of dipolar coupled spin oscillation. To deeply 
understand the coupled dynamics, we calculated the eigenmodes of 
the model system by solving the LLGS equation in a small excitation 
limit. The azimuthal and polar angles of the magnetization in each 
layer have been defined as (pi and respectively. Here / represents 
the layer number of the BPL (/= 1), TPL (/=2), FL (/=3). The LLGS 
equation in Eq. 1 can be described according to the above six 
variables X={xi, ^2, x^, x^, x^) = {(pi, Oi, (p2, O2, ^3, ^3) in the form 

dxi . , , , 

— =Fi[Xi,X2,X^,X4,Xs,Xs) / = 1 ^6 (2) 
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Figure 2 | Resistance and color-coded microwave power spectrum of the 
MTJ cell as a function of the magnetic field, (a) Resistance measured in 
the current density /=7.8X lO^*^ A/m^. Under the negative (positive) field, 
the spins in the free layer are parallel (antiparallel) to the spins in the top 
pinned layer, yielding the lower (higher) magnetoresistance level. The 
spins in the top pinned layer begin to be tilted from the spin flop field 
jUgHsf = 90 mT. The alignments of the spins in the three layers are shown 
as arrow marks for the comparison with the resistance curve, (b), 
(c) Color-coded power spectrums of the MTJ cell in the current biases 
/= + 7.8X10'° A/m" and -7.8X10'° A/m^ respectively, (a) In positive 
current bias, the spins in the free layer are dominantly excited, yielding a 
blue frequency shift with the increasing field, (b) In the negative bias, the 
spins in the top pinned layer are dominantly excited, yielding a red 
frequency shift. A discontinuity in the oscillation mode is observed in both 
cases at around jigH = 45 mT, which is about half the spin-flop field. 



A linearized form of the above equation around the stationary point 
^—{^10^ ^20^ ^30^ •^40' ^50^ ^6o) foflowing Fi{XQ) = 0 can be written as 



dt 



(3) 
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Figure 3 | Microwave power spectrum calculated by on the model MTJ 
device under several different separations (s) between the free layer and 
top pinned layer, (a), (b), (c) Microwave power spectrums as a function of 
the external field calculated from the LLCS equation under the positive 
current bias /= +8.0X10'° A/m" on s = 0.8 nm, 10 nm, 50 nm, 
respectively, (d), (e), (f) The same calculations under the negative current 
bias /= -8.0X 10'° A/m" on s = 0.8 nm, 10 nm, 50 nm, respectively. There 
is also a clear discontinuity and frequency gap in the oscillation mode at 
around 45 mT. As s increases, the gap between the higher frequency mode 
and lower frequency mode decreases. 

Here 6X=X—Xo{3xi=Xi—Xio) represents a smafl deviation from the 
stationary point. If the six eigenvalues and eigenvectors of the 6 X 6 

dFi 

matrix T = Tij = —^ are respectively 8^ and (A;=l~6), i.e., 

TUk = SkUky then 8X in Eq. 3 can be written as 

3X=J2Li^kexp{8kt)Uk (4) 



Here is a numerical constant determined by the initial condition of 
the motion. The motion of each variable is the superposition of the 
six eigenmodes weighted by expie^t). The imaginary part of the 
eigenvalue Im(e^) corresponds to the angular frequency of the 
osciflatory eigenmode, and the real part Re{8k) reveals the decay 
rate of the osciflation. Negative Re{8k) implies a mode whose 
osciflation decays out even though it has an initial osciflation 
amplitude. In the positive Re(£jt), on the other hand, the spin 
osciflation is feasible to be excited^^. In this sense, the value of 
Re(£^) is the crucial factor that determines the damping rate of the 
spin osciflation. 

The matrix f has been numerically calculated under different 
external fields and current biases from the LLGS equation. It gives 
six eigenvalues, three of which are the complex conjugates of the 
others. The three real and imaginary part eigenvalues Re(£^) and 
Im(£^) are depicted as a function of the external field on the positive 
and negative current biases of /= ±3X 10^ A/cm^ in Fig. 4a-4d. The 
mode 1 (A:=l) which has the highest frequency corresponds to 
the excitation of the BPL, and the mode 2 {k=2) and 3 (A; =3) are 
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Figure 4 | Real parts of the eigenvalues of the three coupled magnetic layer system depicted as a function of the external field. (a),(b) Eigenvalues of all 
the three eigenmodes calculated under the positive and negative current biases of /= ± 3 X 10^ A/cm^, respectively. Mode 1 is the excitation of the BPL, and 
the mode 2 and 3 are respectively the optical and acoustic modes between the FL and TPL, which have been identified by the relative phase difference 
between the motions of the three layers, (a) The dominant precession mode where Re(s) is positive changes from the mode 3 (acoustic) to the mode 2 
(optical) at around 45 mT as the field goes higher, (b) The positions of the dominant precession modes are reversed. 



respectively the optical and acoustic modes between the FL and TPL. 
One can note that the optical mode has the higher frequency than the 
acoustic mode because the optical mode has the more dipolar energy. 
Each mode has been identified through the corresponding eigenvec- 
tor as described in the supplementary information A. 

Considering that the spin oscillation is the more heavily excited in 
the mode where Re(£;t) is the larger, the eigenmode with the max- 
imum Re(£;t) among the three modes corresponds to the dominant 
spin precession mode. The dominant precession mode is the acoustic 
(optical) mode at a lower (higher) field of lioH < 30 mT (fioH > 
60 mT) on the case of positive current bias. On the negative current 
bias, it is reversed, i.e., the optical mode is dominant at the lower field 
and the acoustic mode is dominant at the higher field. Compared to 
the coupled motion of the FL and TPL, the motion of BPL is neg- 
ligible. Although BPL is strongly coupled to the TPL by antiferro- 
magnetic exchange coupling, the dipolar field is not fully 
compensated due to the different dynamic motions in the TPL and 
BPL. This reveals that the dipolar interaction should not be neglected 
in the dynamic case. 

The eigenvalue of such a dominant precession mode under several 
different current densities is depicted as a function of the external 
field in Fig. 5. The real parts (imaginary parts) on the positive and 
negative current biases are shown in Fig. 5a and 5b (in Fig. 5c and 5d), 
respectively. The real part representing the decay rate of the oscil- 
lation becomes smaller and goes lower than zero around the field 
45 mT in both current polarities. This agrees with the experimental 
results that the steady state oscillation breaks down around 45 mT. 
The higher current density lifts up the Re(£) values, which can be 
understood as the active reduction of the magnetic damping by the 
current. However, such a general behavior does not occur around the 
breakdown field, where the steady spin precession is inhibited no 
matter how high the current density is. 

The imaginary part depicted in Fig. 5c and 5d shows a frequency 
gap around the breakdown point of 45 mT. The dominant oscilla- 
tion mode changes from the acoustic (optical) mode to the optical 



(acoustic) mode with the increasing field in the positive (negative) 
current bias as depicted in Fig. 5c (Fig. 5d). 

In order to investigate the effect of the dipolar coupling strength, 
we estimated the Re(£) and Im(£) values of the dominant precession 
mode for several different s values. The results of Re(£) as a function 
of external field are shown in Fig. 6a. Figure 6b shows the dependence 
of the frequency gap upon the s values, which was obtained from the 
results of Im(e). The strong reduction of Re(e) around the breakdown 
field observed at small s becomes weaker and weaker as s increases, 
resulting in a negligible dip of Re(£) at s = 1000 nm. The case of s = 
1000 nm can be considered as a system where each layer is coupled to 
the others only through the spin transfer torque without any dipolar 
interactions. In this case, we could not observe any special feature on 
Re(£). This reveals that the anomalous magnetic damping observed 
around the breakdown field is strongly related to the dipolar coup- 
ling. The coupling only through the spin transfer torque cannot make 
such a behavior. 

Recent micromagnetic calculations on a similar spin system 
have considered the dynamic coupling only by the spin transfer 
torque without any dipolar interaction^^'^^. Unlike our results, they 
showed several anomalous features such as the reduction of the 
spectral line-width and frequency jump exactly at the same break- 
down field. But this behavior was observed only at a large current 
density />4.0X 10^ A/cm^. The estimation of eigenvalues and 
eigenmodes assumes a spin motion with small amplitude around 
the stationary points. It cannot cover the dynamics of large ampli- 
tude. In this sense, the anomalous behaviors in purely spin trans- 
fer torque coupled system must come from the nonlinear dyna- 
mics in large amplitude. In fact, we carried out micromagnetic 
calculation at large current density, and it showed the same beha- 
viors as the previous reports as described in the supplementary 
information B. The dipolar- coupling affect the damping property 
even at small spin amplitudes, making the anomalous damping 
behavior be observed experimentally in moderate current range. 
In contrast, the anomalous effect by the dynamic coupling 
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Figure 5 | Real and imaginary parts of the eigenvalue of the coupled spin motion in azimuthal angle on different current densities. The eigenvalues (e) 
and eigenvectors of the coupled spin motion in azimuthal angle have been obtained with the method described in the text. Among the eigenvalues, the 
eigenvalue whose real part is the largest is plotted as a function of external field in the case of s = 0.8 nm. (a),(b) Real part Re(e) calculated on the positive 
and negative current biases, respectively. Re(e) has a tendency to increase with the increasing current densities. However, it becomes smaller and goes 
below zero around 45 mT irrespective of the current densities and polarities. (c),(d) Imaginary part Im(e) calculated on the positive and negative current 
biases, respectively. Im(e) shows similar frequency gap as the experiments and micromagnetic calculations. The relative phase of the spin oscillation 
between the FL and TPL has been identified through the eigenvectors. The schematic diagrams of the corresponding coupled spin motion are depicted in 
the inset of (c) and (d). 



through the mutual spin transfer torque appears only at large spin 
amplitude. 

Damping control through external magnetic field. According to 
the behavior of Re(£) shown in Fig. 5a and 5b, the damping of the 
coupled spin oscillation can be controlled by the external magnetic 
field. The real time motions of the spins in the FL and TPL have been 
obtained from the numerical calculation of Eq. 1 under the following 
condition. Initially the model system was under the field of 27 mT far 
smaller than the breakdown point in the current bias of /= + 
2.0X10^ A/cm^. At the time x = 0 ns, the magnetic field suddenly 



jumped from 27 mT to the breakdown field. The azimuthal angle 
motions of the spins in the FL and TPL have been calculated in 
different s values of 0.8 nm and 20 nm. Figure 7a shows the results 
obtained on the case of s = 0.8 nm. To compare the results with the 
eigenmodes that is valid only small spin oscillation limit, the results 
at a larger x are separately depicted in Fig. 7b. Figure 7c and 7d show 
the results similarly obtained on the case of s = 20 nm. 

As expected from the Re(e) behavior shown in Fig. 5 and Fig. 6, the 
steady state precession below 0 ns decays out above 0 ns. Comparing 
the decaying region (x > 0) with the steady precession region (x < 0), 
one can find that the amplitude ratio between the FL and TPL is not 




Figure 6 | Effect of the dipolar interaction on Re(£) and Im(£) in dominant precession mode, (a) Re(s) of the dominant precession mode depicted as a 
function of the external magnetic field under several different s values. The strong reduction of Re(s) around the breakdown field observed at small s 
becomes weaker and weaker as s increases, resulting in a negligible dip of Re(s) at s = 1000 nm. (b) Frequency gap Af obtained from Im(s) difference 
between the acoustic and optical mode at the breakdown point depicted as a function of s. 
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Figure 7 | Real time motion of the spins calculated in case that the field suddenly changes to the breakdown point, (a), (b) The real time motions of the 
spins in top pinned layer (TPL: black rectangles) and the free layer (FL: red triangles) have been calculated under the field changes described in the text on 
the case of 5 = 0.8 nm. To clarif)^ the beating pattern of the spin oscillation at x>0, dashed blue lines are marked on the points where the ratio of TPL 
amplitude over the FL amplitude has a local maximum. A scale bar representing the inverse of the frequency jump Af^ on the case s = 0.8 nm is drawn to 
be compared with the period of beating pattern, (c), (d) The results obtained similarly on the case of 5 = 20 nm. Af2 represent the frequency jump on the 
case of 5 = 20 nm. 



constant but changes rather periodically. Dashed blue lines marked 
on the points where the ratio of TPL amplitude over the FL amplitude 
has a local maximum show this clearly. This beating pattern becomes 
even clearer when the amplitude of spin oscillation is small enough at 
a larger x, as depicted in Fig. 6b and 6d. The period of beating 
becomes longer as the distance between the FL and TPL becomes 
larger as one can compare Fig. 6d with Fig. 6b. For the comparison of 
the beating frequency with the frequency gap of the eigenmodes, 
scale bars of 1/Afi and 1/Af2 have been drawn in Fig. 7, which 
represent the inverses of the frequency gaps Afj and Af2 obtained 
respectively from the Im(s) values on the two cases s = 0.8 nm and 
20 nm in Fig. 6b. One can find that the beating frequency identified 
by the dashed blue lines coincides with the frequency gap, i.e., the 
frequency difference between optical and acoustic mode at the break- 
down field. This clearly reveals that the beating pattern comes from 
the interference of the two eigenmodes. Considering that the spin 
oscillation in the FL or TPL is expressed as a superposition of the 
acoustic and optical modes by Eq. 4, the observed beating pattern is 
not surprising. The beating pattern appears only near the breakdown 
field because the real part eigenvalues Re(£) of the two competing 
modes are comparable only at this point. Away from this field, the 
difference of Re(s) values is so large that only one eigenmode appears 
dominantly. In this sense, the increased damping at the breakdown 
point is closely related to the beating of the two eigenmodes. 

Since the amplitude of the spin oscillation represents the energy 
stored in the spin motion, the beating pattern of the spin oscillation 
in each layer implies the energy transfer from the FL to TPL, and vice 
versa. In the field away from the breakdown point, both amplitudes 
of the FL and TPL remain the same without any energy transfer 
between the two layers. In the breakdown point, on the other hand, 
the FL and TPL change their amplitudes alternately, exchanging 
energies with each other. 

Recent study on the two dipolar coupled magnetic disks has also 
shown similar energy transfer between the two disks^°. The energy 



transfer in dipolar coupled system is in itself lossless and has no 
relation with the damping of the system. However, if the spin transfer 
torque is involved in the dipolar coupled spin layers like our sample, 
the energy transfer between the coupled layers affects the damping 
property. The current driven active damping on one layer is reversed 
on the other layer. If the FL is excited or anti- damped by the current 
bias, the TPL is the more effectively damped by the same current. 
Thus the transferred energy does not fully come back to the original 
layer due to the intrinsic and current driven active damping of the 
other layer. The energy transferred from one layer with a steady state 
spin precession supported by the current driven active anti-damping 
is more heavily dissipated by the active damping on the other layer. A 
steady spin precession is impossible if only the energy is transferred 
to the other layer. This is the case observed on the breakdown field. In 
spintronic devices like MRAM and STNO, the spin transfer torque is 
considered only to excite the spin state by the active anti-damping. 
But this is true only when the energy exchange between the FL and 
TPL is negligible. At the breakdown field, the FL transfers its energy 
to the TPL through the dipolar coupling, and the current driven spin 
precession ceases to exist. 

Discussion 

To summarize, we observed an anomalous breakdown of the current 
driven spin oscillation in the two spin layers, each of which is anti- 
ferromagnetically coupled to the other by the dipolar interaction, and 
one of which has an additional exchange bias field. From the micro - 
magnetic calculations, the mechanism of the enhanced damping 
could be understood as the energy transfer between the two layers 
and the spin transfer torque acting reversely on each other. This 
damping mechanism gives another way to control the magnetic 
damping. We can control the damping of the coupled system by 
changing the resonant frequencies of the layers, i.e., by changing 
the external field. 
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Besides, the magnetic field acts as a gate to control the energy 
transfer between the two layers. When the field is away from the 
breakdown point, the oscillation energies in the two layers are rather 
isolated. But by putting the field on the breakdown point, the energy 
flows to the other layer. 

Methods 

Sample Preparations. The multilayer film has been grown in the sputtering system 
with the base pressure around 10"^ Pa and post-annealed at 275°C under a magnetic 
field of about 1 T for two hours in order to align the spin orientations of the pinned 
layers in a specific direction. The sample has been patterned into a cylindrical 
nanopillar structure with a diameter of 90 nm by the electron beam lithography 
followed by ion beam etching. Subsequently, the sample was completely covered by 
the thick Si02 film followed by the chemical mechanical polishing (CMP) to open 
only the top contact area. The bottom electrode contact area was patterned through 
the photolithography and ion beam etching up to the bottom layer (TiN). After 
opening the top and bottom electrode contact through an additional 
photolithography, the Ti/Au electrodes have been formed through lift-off process. 
The bottom and top electrode have been formed in a form accessible by GSG type 
radio -frequency (RF) probe. 

Micromagnetic calculations. In the micromagnetic calculation, we did not take into 
account the spin transfer torque between the TPL and BPL. For the calculation of the 
microwave power shown in Fig. 3, the micromagnetic calculation was integrated with 
a model circuit that mimics the realistic measurement circuit. The calculation was 
carried out with the single domain model for each layer, but no significant changes 
were observed from the results using multi-domain model, because the applied 
current density through the MTJ is low enough. The following parameters were used 
fortheMTJ: Layer structure is PtMn layer/ BPL( 1.5 nm)/Ru(0.8 nm)/TPL(1.5 nm)/ 
MgO(0.8 nm~1000 nm)/ FL(2.0 nm). The layer has a circular shape with diameter 
100 nm. The material parameters are as follows. The saturation magnetization = 
1,300 emu/cc, perpendicular interface anisotropy = 1.0 mj/m^ for FL, = 
1240 emu/cc, = 0.9 mj/m^ for TPL, and M, = 1430 emu/cc, exchange bias field 
from the PtMn layer fioHgx =150 mT were used for BPL. Interlayer exchange 
coupling constant between TPL and BPL is — 1.3 |ij/m, damping constant a = 0.01, 
spin polarization P = 0.4, TMR = 30%, and RA value in parallel state of 10 ohm [im^ 
were used. All parameters were chosen to mimic the experiments. In-plane external 
fields were applied along the easy axis with a tilt angle of 10°. Positive external fields 
favor the anti-parallel configuration of magnetizations. For the stochastic calculation, 
the Gaussian-distributed random fluctuation fields (mean = 0, standard deviation = 
[2kBT/(MsVAt)], where At is the integration time step. Vis the volume of unit cell)^^ 
have been added to the effective fields of LLGS equation. 
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